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Abstract 



O 

o 

^ . Using approximately 4.3 million hadronic TP decays collected with the OPAL detector at LEP 

^ '. between 1990 and 1995, we measure the branching fractions of the Z° into up-type and down- 

I type light quarks, and the forward-backward asymmetries, Afb{(1)i using high-momentum 

§^ stable particles as a tag. Adopting a method that employs double tagged events to determine 

' the flavour tagging efficiencies, and assuming the flavour independence of strong interactions 
and SU(2) isospin symmetry, we measure: 



-Rd,s 



0.371 ± 0.016 (stat.) ± 0.016 (syst.) and 



Aj.B(d,s) = 0.068 ±0.035 (stat.) ±0.011 (syst.) , 

assuming the branching fractions and forward-backward asymmetries of down and strange 
quarks to be equal. The results are essentially free of assumptions based on hadronisation 
models. These results are in agreement with the Standard Model expectations and are used to 
infer the left and right handed couplings of strange/down quarks to the Z°, yielding 

g"^^ = -O.Utrol and ^'^ = +0.13t'ofr ■ 

The results for the up quark, Rj{Rd + Rn + Rs) = 0.258 ± 0.031 (stat.) ± 0.032 (syst.) and 
^fb(u) = 0.040 ± 0.067 (stat.) ± 0.028 (syst.), are fully negatively correlated and almost 
completely positively correlated, respectively, with the corresponding down-type results. 
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1 Introduction 



One of the basic assumptions of the Standard Model is flavour universality, namely that apart 
from mass effects, the gauge couplings of all fermions depend only on their charge, weak isospin, 
and colour. While this assumption has been experimentally tested to high accuracy in the 
charged lepton sector, results on quarks are generally less precise. The large number of TP decays 
recorded at the LEP e"'"e~ collider opens the possibility for high-precision measurements of the 
TP couplings to several individual fermion species. Accurate measurements exist particularly 
for the bottom quarks, and to a lesser extent for the charm quark [0. However, at LEP, few 
measurements exist of the couplings of individual up, down, and strange quarks to the neutral 
weak current. The couplings of up and down quarks have been obtained at low from lepton- 
nucleon scattering and from atomic parity violation and are in agreement with the Standard 
Model expectations Direct measurements of light flavours on the Z° resonance are scarce. 
Apart from the overall hadronic width, the yield of photon radiation from quarks is particularly 
sensitive to up- type quarks [|]. A first measurement of the strange quark forward-backward 
asymmetry, which is dependent on hadronisation models, has been published in 0. 

In this paper we determine the decay branching fractions of the Z" into up-type and down- 
type light quarks and the forward-backward asymmetries with only few model assumptions, as 
introduced in 0. The method relies on the property that a high-energy particle in a jet has a 
fiavour correlation with the primary quark. Since charm and bottom quarks do not contribute 
much to the production of long-lived hadrons with a significant fraction of the beam energy, 
a selection of events with ti^, K='=, p(p). Kg, or A (A) baryons of high energy provides rather 
pure samples of light fiavoursQ. The main challenge is to determine the relative yields of up, 
down, and strange quarks in such event samples. This is achieved by using the information 
from double tagged events, where two high-momentum particles are found in opposite event 
hemispheres, together with some general hadronisation symmetries. 

The method of measuring the branching fractions and the forward-backward asymmetries is 
described in detail in Section 2. The elements of the OPAL detector pertinent to this analysis 
and the basic event selection are presented in Sections 3 and 4. In Section 5 we describe 
the estimation of the charm and bottom backgrounds in our high-momentum samples and in 
Section 6 the electroweak observables are determined for up-type and down-type light quarks. 
Systematic uncertainties are evaluated in Section 7 and the conclusions are given in Section 8. 



2 The Method 

The observed long range charge correlations in e^e~ events |^ show that the particle with the 
highest energy in a jet tends to carry the quantum numbers of the primary quark. Thus, jets 
from primary strange quarks lead to high-energy strange hadrons, for example, and primary up 
and down quarks lead to high-energy pions and protons. We base our analysis on this property 
and tag light quark events by identifying high-energy charged pions, charged and neutral kaons, 
protons, and A baryons with Xp = 2ph/Ecm > 0.5, where ph is the momentum of the tagging 

^Unless explicitly stated otherwise, charge conjugation of the tagging particle and primary quark types is 
implied throughout this paper. 
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hadron, h, and E^m is the centre-of-mass energy of the event. The value of Xp > 0.5 is chosen 
to minimise the statistical and systematic uncertainties. We now make a general discussion of 
the method without considering the detector. 



2.1 Method to Determine the Branching Fractions 

If the relations between particle types and primary flavours were unambiguous, double tagged 
events, in which a tagging hadron is found in each of the event hemispheres, could be used to 
determine the flavour tagging efficiencies in a straightforward way. The method would then be 
essentially free from uncertainties due to the detailed properties of the flavour tag and detector 
effects. Almost no ambiguity exists for bottom particles, where the double tagging method has 
been successfully applied (see [§], for example). On the other hand, high-energy light-flavour 
mesons and baryons can be produced by processes other than the hadronisation of primary 
quarks. Leading mesons contain a quark and an antiquark, either of which could be the quark 
into which the Z° decayed directly, which introduces some ambiguity. In addition, decays of 
these so-called primary hadrons and of particles produced from the hadronisation sea tend to 
obscure further the primary flavour source of each particle type. Thus, each particle species is 
produced from a mixture of several primary quark flavours. Neglecting correlations, the number 
of tagged event hemispheres (as defined by the thrust axis) and the number of double tagged 
events can be expressed as: 



N, 



h 



had 



2 Yl V'R, and (1) 

g=d,u,s,c,b 



= E (2) 

^Vhad g=d,u,s,c,b 

where Nh is the number of hemispheres with a tagging hadron h = vr^, K^, p(p). Kg, A(A), 
A'had is the number of hadronic Z° decays, and N^h' is the number of double tagged events with 
tagging particle types h and h'. The rj^ denote the fraction of hemispheres with a primary 
quark fiavour q which are tagged by a hadron of type h, and Rg is the hadronic branching 
fraction of the Z° to quarks q: 

R, = (3) 

J- had 

As will be shown in Section 5, the charm and bottom fractions can be determined separately 
in a straightforward way, leaving fifteen unknown rj^ and three unknown Rq. On the other hand, 
for the five tagging particle types we consider, there is a system of five equations for tagged 
hemispheres (Eq. 1) and fifteen equations for double tagged event types (Eq. 2). Due to the non- 
linearity and degeneracy of the equation system, it is not solvable and additional constraints 
have to be found to obtain a solution. 

As discussed in detail in 0, the necessary constraints can be derived from SU(2) isospin 
symmetry and the fiavour independence of QCD. We use the following hadronisation relations: 

'if = <* , (4) 
= (5) 
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vT^y = Vf, (6) 
r/s^"^^"^ = r,r,and (7) 
= (8) 



Note that K°(K°) implies Kg plus Kl- Small deviations from these relations are discussed later 
in Section 7. In addition, the overall normalisation J2Rq = ^ provides another constraint: 

+ Rd + Rs = I - Rc - Rh = 0.620 ± 0.010 , (9) 

where the LEP measurements [0] of i?b and Rc can be used to constrain the sum of the light- 
flavour branching fractions. To be independent of the measurements of the heavy flavour 
fractions, we also express our results in terms of: 

r^jd + + Fss Rd + Ru + Rs 

The R' are related to the R„ via Rq — Rqi^ — Rh — Rc) ■ 

These constraints are still not sufficient to allow the system to be solved. Motivated by 
the weak isospin structure of the Standard Model, we assume Rd = Rs, thereby reducing the 
number of unknown Rg and making possible a solution of the equation system. 



2.2 Method to Determine the Forward-Backward Asymmetries 

The other directly measurable electroweak-related observable is the forward-backward 
asymmetry Ap^ of a hadron h: 

,h _ iV^(cosg>0) - iVfe(cosg<0) 

^FB — Tt 5 {^^) 

where 9 is the angle of the tagging hadron h with respect to the incoming electron direction. 
The relation between this directly observable asymmetry and the desired forward-backward 
asymmetries of the quarks is given by 0: 

= E {^./.'(2rJ - l)}AMq) , (12) 
q 

where the reliabilities of the charge tag are given by = N^""^^^ /Nh^g, iV^°j''S'=t being the 
number of hadrons h which have the same sign of the charge as the primary quark q. The 
reliability therefore takes into account dilutions due to the misidentification of the sign of the 
charges of the quarks, Sg. Finally, fg is the fraction of tagged hadrons h stemming from a 
primary quark q, i.e. = (rj^Rg)/ (J2q' Vq'^g')- Therefore, the forward-backward asymmetries 
Apsiq) can be determined only after the r]^, Rg, and reliabilities are known. The reliabilities 
can be determined from the ratio of double tagged events with tagging particles of opposite 
charge, N^/^^ over the total number of double tagged events, Nhh'-. 

= E 7yM^. W + (1 - --JXl - <■)} . (13) 

^^hh' g [l^q' Vq'Vq' -ttq' ) 
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Since the Kg provides no information on the charge of the primary quark, we restrict the 
determination of the forward-backward asymmetries to charged pions, charged kaons, protons, 
and A baryons. With these four hadron types we obtain four measurements of Ap^ and ten 
ratios N^y^ /Nhw which must be used to determine 12 unknown hght flavour reliabilities and 
the three asymmetries ApB{q)- Here we assume that the heavy flavour terms can be determined 
separately, as will be shown later in Section 5. As for the solution of the equation system for 
the branching fractions, we invoke hadronisation symmetries based on SU(2) isospin invariance 
to reduce the number of unknown reliabilities: 



rl = rl , (14) 
rf = 0.5 , and (15) 
r^(^) = l-r^^'^K (16) 

Note that for the baryons (the proton and A), the tagging hadron carries the sign of the baryon 
number of the primary quark, not the electric charge, so that r^'^''^"'^ < 0.5 for down-type quarks 
and r^^''y°'^ > 0.5 for up-type quarks. In addition, in order to solve the equation system we 
must assume that (see 0): 

r^* = 0.20 ±0.10 , (17) 

where the value is taken from the JETSET model and is assigned a large error to take into 
account uncertainties in the JETSET modelling. This particular reliability is chosen as the one 
to be fixed since is small compared to the other flavour fractions. Therefore, even a large 
uncertainty in does not affect significantly the final result. 

In the preceding discussion we have ignored biases due to geometrical and kinematic 
constraints which have to be taken into account. Requiring a high-energy hadron in an event 
reduces the phase space for gluon bremsstrahlung and thus introduces a kinematic correlation 
between the event hemispheres. Similarly, restrictions on the geometrical acceptance introduce 
corrections to the equation systems we use. We parametrise this correction by a factor p, so 
that the double tagging probability, rj^^ , is given by: 

= Pvivi' ■ (18) 

In the ideal case of no bias, p = 1. Model calculations suggest that this correlation is 
independent of the tagging particle type and the primary quark flavour and is approximately 
1.07 for a sample of tagging hadrons with Xp > 0.5 without experimental cuts and neglecting 
detector effects, and is essentially due to gluon bremsstrahlung. In order for the measurement 
of the light-flavour electroweak parameters to be as model independent as possible, we let p 
be a free parameter in the equation system but assume that it is independent of the tagging 
hadron species. Note that p is largely uncorrelated with the other parameters. 



3 The OPAL Detector 



The OPAL detector is described in detail in [0. Here we summarise only the features of the 
detector which are important for this analysis. 
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Central to this study is the determination of the momentum and the identity of different 
types of charged hadrons and neutral hadrons which decay to stable charged particles. OPAL 
has a system of tracking devices inside a solenoid which provides a magnetic field of 0.435 T. A 
charged track momentum resolution of (Jp/p = 0.02©0.0015 pt, where pt is the component of the 
total track momentum p in the plane perpendicular to the beam axis measured in GeV, has been 



achieved. The innermost part is a silicon microvertex detector |11], surrounded by three drift 
chambers: a vertex detector, a large volume jet chamber which provides up to 159 space points 
per track, and z-chambers to give a more precise measurement of the polar angle of charged 
track^. The large number of measurements in the jet chamber also provides a determination of 



the specific ionisation energy loss, dE/dx, with a resolution of a{dE/dx)/{dE/dx) ~ 0.035 [|T2 
for well-separated tracks with | cos6\ < 0.7. This resolution allows the identification of charged 
pions, charged kaons, and protons up to the highest particle momenta ||13[. The large radius of 
the jet chamber {R = 185 cm) also allows a high reconstruction efficiency for large-momentum, 
weakly decaying hadrons with relatively long decay lengths, such as Kg — > 7r"'"7r~ and A —>■ pvr". 

To estimate the contributions from bottom quarks, we also identify secondary vertices from 
b hadron decays |T^ and use electron and muon identification |T^, to tag semi-leptonic 



bottom decays. The identification of a secondary vertex profits particularly from the high- 
precision measurements in the silicon microvertex detector. The lepton identification is largely 
based on the electromagnetic calorimeter, which consists of 11704 lead glass blocks which 
subtend a solid angle of 40 x 40 mrad, and muon chambers which are placed behind an average 
of eight absorption lengths of detector material. 

To determine detector efficiencies and possible detector biases, we use a sample of 
approximately eight million simulated hadronic Z° decays generated with the JETSET model H] 



and passed through a detailed simulation of the OPAL detector jl^. The fragmentation 



parameters have been tuned to describe overall event shapes and distributions as described 



in detail in |]T8|. Events are generated with two versions of JETSET. Events with JETSET 7.3 
use a simulation of the detector up to and including 1993, while the JETSET 7.4 events use a 
simulation of the detector corresponding to subsequent years. 



4 Event Selection 



The analysis is based on approximately 4.3 million multihadronic Z° decays collected between 
1990 and 1995. Of these events approximately 90% were collected at centre-of-mass energies 
within ±200 MeV of the Z° mass, and the rest of the events within ±3 GeV above and below 
the Z° peak. The standard OPAL requirements for the multihadronic event selection are given 
in To enrich the light-fiavour fraction we retain only those events which have at least 

one well measured charged track, a Kg, or a A baryon, with a scaled momentum Xp > 0.5. 
Details of the particle selections are given below. Since this high-a;^ cut rejects relatively few 
Z° — >■ r+r" events in our sample, we further demand at least eight good charged tracks 



in an event. In addition, to ensure good vr^, K^ and proton separation, and reliable Kg 

^ OPAL uses a right handed coordinate system in which the z axis points along the direction of the electron 
beam, r is the coordinate normal to this axis, and 6 and (j) are the polar and azimuthal angles with respect to 
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and A reconstruction, we require that the tagging particles have a polar angle | cos^| < 0.7. 
For |cos^^| > 0.7, the dE/dx separation of kaons and protons from pions degrades rapidly 
with increasing | cos^^| of the tracks. In order to select events which are well contained in the 
detector, we restrict the polar angle of the thrust axis (calculated using both charged tracks and 
clusters in the electromagnetic calorimeter which have no associated track in the jet chamber) to 
satisfy | cos ^Thrust I < 0.8. With these requirements, 198 309 events are retained. The remaining 
background from Z° t^t~ events in our sample is negligible (less than 0.05%), as estimated 
from Monte Carlo simulations. 

Crucial to this analysis is the identification and separation of samples of charged pions, 
charged kaons and protons, which is achieved using the dE/dx measurement, and Kg and A 
baryons, which are observed by reconstructing their decay vertices and calculating the invariant 
mass of the decay products. Efficiency losses which are common to all particle types, such as the 
finite geometric acceptance, etc., can be absorbed into the values of 77^. However, other sources, 
such as those due to particle identification {dE/dx, secondary vertex finding) have to be taken 
into account for the hadronisation symmetries. Similarly, there exists some misidentification 
probability leading to a migration from the true to an apparent particle identity, which is 
especially relevant for the stable charged particles. This must also be taken into account and 
the hadronisation symmetries corrected for different detection efficiencies and migrations. Due 
to the fundamentally different reconstruction of the weakly decaying Kg and A, these samples 
are largely decoupled from the charged hadron samples. The relation between apparent particle 
type /idet and true particle type h is given by some fiow matrix, Ej^^^, such that 

Nh.. = Y.^LNh . (19) 

h 

Similarly, the forward-backward asymmetries for the measured samples A!^^ are related to the 
pure particle type asymmetries A^^ by 

^f^'b* = T.^LAb ■ (20) 

h 

Apart from reducing the discrimination power between the various fiavours, some systematic 
uncertainties related both to the efficiency and purity are introduced. 

For the dE/dx measurement, tracks are required to have at least 20 hits in the jet chamber 
used in the calculation of the energy loss, a polar angle 9 satisfying | cos^| < 0.7, a distance 
of closest approach to the interaction point in the plane orthogonal to the beam direction of 
|(io| < 2 mm, and the corresponding distance along the beam direction \zq\ < 40 cm. To 
suppress badly measured tracks we further restrict the scaled momentum to Xp < 1.07, which 
takes into account the 7% momentum resolution of a track with Xp = 1. A study of muon pair 
events finds that less than 0.3% of tracks which have the full beam energy have an incorrect 
charge assignment. 

To separate samples of pions, kaons and protons, we use for each track the dE/dx 
weight, Wh, which is defined as the signed probability of the track to be consistent 
with a certain particle species hypothesis, h. The sign represents the sign of the difference 
between the measured d£^/dx(meas.) and the expected di?/da;(exp.) for the particle species 
hypothesis, namely a positive weight if d-E/da:(meas.)>d£'/dx(exp.) and a negative weight if 
d£'/da;(meas.)<d£'/da;(exp.). In particular, we require: 
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• for pion candidates: {Wt,± > 0.01 or < —0.1), and \wk± \ < 0.1 ; 

• for kaon candidates: \wk± \ > 0.1 and |w7r±| < 0.1 ; 

• for proton candidates: (wp(p) > 0.1 or Wp(p) < —0.01), and \wk± \ < 0.1 . 

These selection criteria result in three disjoint track samples. The dE/dx separation power as a 
function of charged track momentum is shown in Fig. 1. In addition, charged tracks are rejected 



if they pass either the standard electron or muon identification requirements 1T^,|T6| . After these 
requirements, backgrounds from electrons and muons are negligible, as determined from Monte 
Carlo simulations. The largest contamination is the 0.8% of pion candidates which are actually 
muons. In addition, some contamination from charged hyperons (mostly S~) is present in 
the proton sample at the 10% level according to the Monte Carlo simulations. However, since 
we make no hadronisation assumptions about the proton sample, this contamination is not a 
problem for the analysis. 

These cuts (including the event and thrust axis cuts) lead to efficiencies of (34.7±0.1±5.2)% 
for pion tagged hemispheres, (29.5 ± 0.1 ± 4.6)% for kaon tagged hemispheres and (23.4 ± 0.3 ± 
5.3)% for proton tagged hemispheres, where in each case the uncertainties are statistical and 
systematic. The efficiency is defined as the number of event hemispheres which are tagged at the 
detector level divided by the number of event hemispheres which are tagged at the Monte Carlo 
generator level. Differences in the tagging efficiencies are due entirely to the dE/dx selections. 
The systematic uncertainty is determined by scaling the corresponding widths of the dE/dx 
distributions, a, by 1.00 ±0.05 and by varying the expected dE/dx value for each particle type 
by ±0.15(7. The ranges of these uncertainties come from checks with pions from Kg tt+tt" 
decays and protons from A pvr^ decays for track momenta greater than 2 GeV in which 
the means and widths of the normalised dE/dx distributions showed maximum deviations of 
±0.15(7 and ±5%, respectively. The uncertainties are therefore strongly positively correlated 
between the various particle species. The other selection criteria were studied in the Monte 
Carlo event samples and found to have no significant bias to select any one tagging particle 
type over another. 

The procedures to identify the weakly decaying hadrons Kg and A are described in detail 
in ||21[] and ||22|, respectively. Here we summarise the main ingredients. To find Kg —>■ tt+tt" 



and A — > pvr" candidates, we combine two oppositely charged tracks which have at least 20 hits 
in the jet chamber. We then search for a crossing point of these tracks in the plane orthogonal 
to the beam axis. If a good secondary vertex is found, the vr+vr" and pvr" invariant masses of 
the combinations are calculated. Good Kg candidates are required to have Xp > 0.5 and have 
invariant masses in the ranges 430 MeV < m^+^- < 570 MeV and rrip^- > 1.13 GeV, in order to 
reduce the contamination from A —>■ pvr^ decays, where is the invariant mass of the tracks 

assuming that they are a pair of pions and mp^^- a proton-pion pair. Likewise, all candidates 
which have Xp > 0.5 and 1.10 GeV < mp^- < 1.13 GeV are accepted as high-momentum 
A candidates. No additional rejection of Kg —>■ tt+tt^ background in the A sample using an 
invariant mass cut is made since such a rejection cannot be done efficiently. The invariant 
mass distributions for all candidates which pass these selection criteria are shown in Figs. 2 
and 3. The detection efficiencies are determined using the Monte Carlo simulated events. After 
correcting for mass resolution differences between data and Monte Carlo simulation as in ||22 



the efficiencies, defined as before, are found to be (9.4 ± 0.1 ± 0.3)% for Kg tagged hemispheres 
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and (4.7±0.1 ±0.3)% for A tagged hemispheres, where the errors are statistical and systematic, 
respectively. The systematic errors are mainly due to uncertainties in the correction factors 
for the different mass resolutions in the data and Monte Carlo and the simulation of other cut 
distributions. Detailed treatments of the determination of these systematic errors can be found 
in |21|,|22|. With these requirements we select the number of tagged hemispheres and double 
tagged events listed in Table 1. 



The purities of the final samples, as determined from the efficiencies calculated above 
with appropriate weighting by the relative particle yields in the data are given in Table 2. 
Although there is some uncertainty in the particle identification efficiencies, the charged pion, 
kaon and proton sample purities are known to ±0.011, ±0.013 and ±0.052, respectively. The 
errors are smaller for the purities than for the efficiencies because of the positive correlations 
among the efficiencies. The Kg and A sample purities are largely uncorrelated with those of 
the charged hadron samples. In addition to small combinatorial backgrounds composed of 
randomly paired tracks. Kg tt+tt^ decays form the principal background in the A sample, 
and A —>■ pyr^ combinations form a somewhat less important background in the Kg sample. The 
total backgrounds are determined by fitting to the invariant mass spectrum shown in Fig. 2, as 
in ^T|, and by using a sideband method to determine the background in the distribution shown 
in Fig. 3, as in ||22|. For the Kg sample, after subtracting the contribution from A — pvr" decays. 



a combinatorial background level of (5.6 ± 2.9)% is found, where the error is systematic ET 



Combining this with the relative uncertainty on the efficiency, the total relative systematic 
error on the Kg yield is ±4.1%. For the A, the combinatorial background present in the sample 
is (7.7 ± 2.4)%, where the error is again systematic p2[. The total relative error on the A yield 
is therefore ±6.3%. 



hadron 


tagged 


double tagged events 


type 


hemispheres 




K± 


P(p) 


J^s 


A(A) 




52170 


392 


416 


99 


46 


18 


K± 


40 229 




265 


136 


48 


15 


P(P) 


9 350 






13 


15 


4 


K°s 


5 026 








3 


1 


A(A) 


1349 










1 



Table 1: Number of tagged event hemispheres and double tagged events for Xp > 0.5. 



5 Heavy Quark Contributions 

Before determining the light-flavour electroweak properties, we flrst flnd the contributions 
of charm and bottom quarks to the selected samples. In the case of the bottom quark, we 
determine the fraction directly from the data (except for the tagged A sample) by searching 
for a secondary vertex or for a lepton with a high transverse momentum relative to the jet 
direction in the high-Xp samples. Lack of statistics prevents a similar method from being used 
to determine the charm backgrounds with the necessary precision. Instead we use the JETSET 
model and estimate the uncertainties based on the measured production and decay properties 
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assigned 


true TT^ 


true 


true p(p) 


true Kg 


true A (A) 


7r± 


0.895 


0.080 


0.002 


0.011 


0.000 


K± 


0.106 


0.712 


0.171 


0.001 


0.008 


P(p) 


0.013 


0.367 


0.591 


0.001 


0.028 


Kg 


0.003 


0.001 


0.000 


0.979 


0.017 


A(A) 


0.001 


0.003 


0.000 


0.307 


0.690 



Table 2: The particle composition in samples assigned as charged pions, charged and neutral 
kaons, protons, A baryons. Note that in some of the rows the numbers do not add up to 100%, 
due to the presence of small backgrounds from other sources such as leptons. The combinatorial 
backgrounds under the Kg and A signals are not considered here and are discussed in the text. 
Charged hyperons were treated as "protons" for the purposes of identification. 



of charmed hadrons. Due to the low statistics available, the b quark flavour fractions and 
reliabilities for the tagged A sample are also taken from the simulation. 



5.1 Bottom Quark Contributions 

To determine the fractions //^ of Z" — > bb events in the high-Xp event samples, we search for a 
secondary vertex displaced from the primary event vertex in the event hemisphere opposite to 



that of the high-Xp tagging particle. The details of the vertex finding are discussed in \T^\. We 
require the vertex to have at least four assigned tracks and to have a decay length significance 
l/ai >10, where I is the observed distance to the event vertex and ai its error. This yields a b 
tagging efficiency per event hemisphere, eb, of about 16% with a purity, Pb, of (96 ± 1)%, using 
the methods in The desired fractions are then given by: 



and are listed in Table 3. Here A^^~*^^ is the observed number of events with a tagged hadron h 
in one hemisphere and a secondary b vertex in the one opposite and A^^~*'^§ is the total number 
of hemispheres with a secondary b vertex. The term Bh takes into account the contribution to 
the measured fraction from events tagged as a b event but which do not come from Z° — bb 
events. This term, which is weighted by (1 — Pb), is calculated from Monte Carlo simulated 
events and is the only model-dependent number used in the measurement of the backgrounds 
from Z° — > bb events. A conservative error of 50% is assigned to the calculated value of Bh, 
which for pions, for example, is calculated to be i3^± = 0.127 ± 0.015 ± 0.064, where the 
first error is statistical and the second systematic. The values for the other hadron types are 
similar. Given the very high purity of the b tagged sample, the infiuence of this term on the final 
measurement is small. However, other sources of systematic uncertainty, such as uncertainties 

in Pb and Pb, are negligible in comparison. Due to low statistics, f^^^^ is taken from the 
simulation with a ±50% relative systematic error, as given in Table 3. 

In order to verify that the method does not introduce systematic biases, the measurement 
is repeated on JETSET events with a detailed simulation of the OPAL detector and the results 
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are compared with the true b fractions in the Monte Carlo event sample. No systematic biases 
are found. The measurement is also checked with the data by selecting b events using leptons 
with a high transverse momentum with respect to the jet direction (described in more detail 
below). The two measurements are found to be in good agreement, although the statistical 
errors for the lepton samples are large. The dependence of the measurement on the purity of 
the b tagged sample is also checked by varying the minimum decay length significance over the 
range from 6 to 14. No systematic dependence of the measured b fractions is found. 



hadron type 


ft 




7r± 


0.078±0.004 


0.79±0.11 


K± 


0.039±0.004 


0.67±0.12 


P(P) 


0.051±0.009 


0.71±0.30 


Kg 


0.036±0.010 




A(A) 


0.031±0.015 


0.29±0.11 



Table 3: Fractions and reliabilities of hadrons from bottom events, where the errors include 
both statistical and systematic uncertainties added in quadrature. 



To obtain the high-Xp hadron charge reliabilities in bottom events one has to know the charge 
of the tagged bottom hadron. Instead of secondary vertices, we therefore use high-energy muons 
to identify bottom particles. Details of the muon identification and tagging procedures are given 
in [T^. The observed charged tracks and neutral clusters in the electromagnetic calorimeter 
are combined into jets of a maximum observable invariant mass of 7 GeV using the JADE 



algorithm with the 'EO' recombination scheme ||2^. The event is tagged as a Z — > bb event 



if a muon candidate is found with a total momentum greater than 3 GeV and a transverse 
momentum with respect to the jet axis, px, of at least 1 GeV. 

First we determine the overall charge reliability of the high-pT muon using double tagged b 
events. The b reliability is given by: 




(2rb-l) = W2-^-l, (22) 



where A^bb is the total number of double tagged b events and A^btf ^ is the number of such events 
in which the tagging leptons have opposite charges. We measure rb = 0.650 ± 0.015, where 
the error is statistical. Note that this measurement includes effects such as B^B mixing. The 
high-Xp hadron charge reliability is then given by: 

{2r^ - l) = , (23) 

where j\r^~*'^S'*-'^P jg the number of events in A^^~*^^ in which the tagging particles have opposite 

charges. The results are also listed in Table 3. Again, due to low statistics, r^^^^ is taken from 

the simulation with a systematic error corresponding to a relative error of ±50% on (2r^^^^-l), 
which is the factor relevant for the measurement of the forward-backward asymmetries. 
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5.2 Charm Quark Contributions 



In the absence of an efficient enough charm tag we estimate the fl^ from Monte Carlo 
simulations. The estimated uncertainties on the /^'^ are based on the knowledge of charm 
production and decay. The essential ingredients are the fragmentation functions of charm 
quarks and the relative production yields of charmed hadrons ||2^ and their decay 
properties 0. Taking into account the corresponding measurements and their uncertainties, 
we derive the fractions as listed in Table 4. 



To cross-check these estimates with the data, we search in the tagged samples of high-Xp 
particles for D*^ mesons in the opposite hemisphere. We use the transition D*^ —>■ D^vr^ and 
the D° decay modes D° K"7r+, K~e"'"z/e, and K~yU+z/^. Details of the selection procedures are 
given in [^]. The results are consistent with the above estimates although with large statistical 
uncertainties. Results are reported for comparison in Table 4. 



hadron type 


f!^ (JETSET) 


/c' (DATA) 


(JETSET) 




0.068±0.007 


0.04 ±0.04 


0.70±0.12 


K± 


O.lOliO.Oll 


0.17 ±0.06 


0.25±0.06 


P(p) 


0.088±0.009 




0.40±0.12 


Kg 


0.114±0.012 






A(A) 


0.091±0.010 




0.70±0.12 



Table 4: Fractions and reliabilities for hadrons from charm events, where the errors include both 
statistical and uncorrelated systematic uncertainties added in quadrature. No measurements 

were possible for /^®, /cf^, and f^^^^ in the data due to lack of statistics. In addition to these 
uncertainties in the individual particle fractions, overall shifts between —5% and ±20% are 
observed due to uncertainties in the charm fragmentation function and have to be taken into 
account. 



We also take the high-Xp hadron charge reliabilities in charm events from studies of simulated 
events. The values are also given in Table 4. As a cross-check we measure the charge reliabilities 
from data using the D* sample described above, using Eq. |23|, but taking the overall charm 
reliability Vc from simulation. The results are compatible with the values from the simulation 
but have large statistical errors. 

Several sources of systematic uncertainties have been considered. Of particular relevance are 
those affecting the relative fractions of Kg and from charm decays. The ratio of Kg to K^ 
production depends on the inclusive branching ratios of the D°, D^*", Dg and A+ into the tagged 
particle species and the corresponding decay multiplicities. The former were taken from 0] 
and propagated through to the fl^. The decay multiplicities as used in the JETSET simulation 



were compared to the measurements of the Mark III collaboration p8[. The differences lead to 
uncertainties of ~4% for each particle species. The resulting errors for each of them are given 
in Table 4. In addition, the of all hadron species are affected coherently by the hardness 
of the fragmentation function. Allowing the Cc parameter of the Peterson et al. fragmentation 
function |29| to vary between 0.020< ec <0.039, which corresponds to 0.474 < (xg) < 0.501, 
as given by the measurements of the D* and D fragmentation functions overall changes of 
—5% to ±20% are introduced. 
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6 The Branching Fractions and Forward- Backward 
Asymmetries of Light Flavours 



As discussed in Section 2 the branching fractions of the Z° into hght quarks and the forward- 
backward asymmetries are determined under the assumption that the electroweak couphngs for 
down and strange quarks are the same. We start with the branching fractions and afterwards 
discuss the forward-backward asymmetries in Section 6.2, considering at first only statistical 
uncertainties. The systematic uncertainties are addressed in Section 7. 

6.1 The Branching Fractions R'^ 

After subtracting the contributions from heavy quarks, we solve the equation system for -Rd,s- 
In addition to the five measured numbers of tagged hemispheres and fifteen combinations 
of double tagged events (Eqs. 1—2), we use the five hadronisation symmetries (Eqs. 4—8). 
Therefore, there are twenty measurements and twelve unknown parameters, including -Rd,s, a 
global p, and ten 77^. The equation system is solved by minimising a function which uses 
the 20 measured quantities of Eqs. 1—2. In order to minimise the dependence of our results 
on heavy flavour electroweak measurements, we determine the partial light flavour branching 
fractions: 

i?d,s = 0.371 ±0.016, (24) 

where the error is purely statistical and the value of is 4.4 for eight degrees of freedom at the 
minimum. The result for R'^ = 0.258 ± 0.031 is fully anticorrelated to the value of R'^^ since 
R'^ + R'^ + R'^ = 1. The solution is tested for stability and uniqueness by using different starting 
values for the unknown parameters. No other physical solution is found. The Standard Model 
expectations for a top quark mass of 175 GeV and a Higgs boson mass of 300 GeV are 0.359 
and 0.282, respectively, consistent with the above measurements. If we use the world average 
values for Rc and i?b 0? we can determine the branching fraction of the Z° into down-type 
light quarks: 

i?d,s = 0.230 ± 0.010 , (25) 

where again the error is purely statistical. The corresponding _Ru = 0.160 ± 0.019 is fully 
anticorrelated with this result. 

The values for the other parameters are given in Table 5, along with those obtained from 
the solution using Monte Carlo events for comparison. Some differences are found for r/J and 
for the baryons. The value of the hemisphere correlation, p, is comparable in data and Monte 
Carlo. Most of the correlation is due to the limited geometric acceptance and the selection 
criteria, which are well modelled in the Monte Carlo simulations. The correlation in JETSET 
events without detector simulation and without any restrictions on the geometric acceptance 
is about 1.07 and is due to hemisphere correlations from gluon radiation |Q. 
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parameter 


OPAL 


JETSET (input) 


P 


1.290 ± 0.043 ± 0.025 


1.286 ± 0.032 (1.238) 


vf 


0.0406 ± 0.0016 ± 0.0021 


0.0441 ± 0.0012 (0.0438) 


vf 


0.0074 ± 0.0028 ± 0.0011 


0.0084 ± 0.0021 (0.0089) 




0.0010 ± 0.0052 ± 0.0012 


0.0038 ± 0.0024 (0.0045) 




0.0230 ± 0.0052 ± 0.0011 


0.0178 ± 0.0025 (0.0112) 




0.0344 ± 0.0029 ± 0.0013 


0.0363 ± 0.0027 (0.0407) 




0.0006 ± 0.0027 ± 0.0008 


0.0065 ± 0.0032 (0.0042) 


r)P(p) 
/u 


0.0082 ± 0.0056 ± 0.0015 


0.0086 ± 0.0048 (0.0124) 


^p{p) 


0.0086 ± 0.0023 ± 0.0005 


0.0042 ± 0.0017 (0.0037) 




0.0049 ± 0.0020 ± 0.0004 


0.0032 ± 0.0014 (0.0022) 


^A(A) 


0.0040 ± 0.0033 ± 0.0005 


0.0107 ± 0.0024 (0.0129) 



Table 5: Free parameters in the equation system and their fitted values, with statistical and 
systematic errors, which are highly correlated. Also given for comparison are the values from 
a solution using JETSET events with full detector simulation with statistical errors only and 
the input values given in parentheses. Good agreement between the Monte Carlo input values 
and solutions is found. 

6.2 The Forward- Backward Asymmetries 

To determine the forward-backward asymmetries, we divide the data sample into events 
collected within ±200 MeV of the Z° peak with an average centre-of-mass energy {Ecm) = 
91.2 GeV, and those below and above the peak, which have average centre-of-mass energies 
of (-Ecm) = 89.5 GeV and (-Ecm) = 92.9 GeV, respectively. The backgrounds due to charm 
and bottom events are known since the heavy flavour terms in Eq. 12 have been measured 
separately^]. We then use Eq. 12 together with the relation of Eq. 13 to determine the reliabilities 
of the charge tags for the light flavours. The measured hadron forward-backward asymmetries, 
which have been determined by maximising a log-likelihood function, are listed in Table 6. 



hadron type 






Ecm > 


7r± 


(-1.22 ±2.14) X 10"^ 


(-1.60 ±0.52) X 10-2 


(±0.12 ± 1.78) X 10^2 


K± 


(±1.76 ±2.52) X 10-2 


(-3.12 ±0.58) X 10-2 


(-6.66 ±2.04) X 10-2 


P(p) 


(±9.98 ±5.18) X 10-2 


(-1.28 ± 1.26) X 10-2 


(±3.42 ±4.42) X 10-2 


A(A) 


(±10.0 ± 12.8) X 10-2 


(±4.32 ±3.18) X 10-2 


(±8.8 ± 11.4) X 10-2 



Table 6: The measured forward-backward asymmetries, ^4^^*^', with statistical errors for various 
centre-of-mass energies, Ecm- Note that the A^§^ have not been corrected for the purities of 
identified samples. Results are given in percent. 



We extend the assumption = Rs to the equality of the forward-backward asymmetries 
of down and strange quarks, and use the hadronisation symmetries of Eqs. 14—16 and the 

■^See Tables 3-4 for the heavy flavour fractions and rehabiUties r^, and |^ for Apb{c) and ApB{h). 
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constraint from Eq. 17. This leaves 8 unknown reliabilities, 10 combinations of double tagged 
events of same and opposite charges (note that double tagged events with cannot be used), 
and 4 measured hadron A^^pj^ with which to determine the light flavour forward-backward 
asymmetries. We find for energies near the Z° peak: 

^FB(d,s) = 0.068 ±0.035 and Afb{v) = 0.040 ±0.067, (26) 

where the errors given are statistical only. The correlation between these results is +91%. 
The rehabilities are found to be consistent with the JETSET predictions, although with large 
statistical errors. The three best measured reliabilities are rf = 0.81±0.02, = 0.70±0.18, 
and = 1.05 ± 0.06, where the errors are statistical. Near the Z'^ peak, the flavour 

composition changes very slowly and QCD scaling violations are likewise negligible. This is 
consistent with results from studies using JETSET where the flavour fractions and reliabilities 
for the off-peak data are found to be the same as for the on-peak data. Therefore, we use 
the flavour fractions and reliabilities as determined from the on-peak data to determine the 
off-peak forward-backward asymmetries. We find for energies above the Z'^ peak: 

AFB{d,s) = 0.228 ±0.082 and Afb{vl) = 0.397 ±0.183, (27) 

and below the Z° peak: 

AFB{d,s) = 0.019 ±0.098, and Afb{vl) = 0.103 ±0.216. (28) 

The values as a function of the centre-of-mass energy are shown in Fig. 4, along with the 
predictions of the Standard Model, which are found to be consistent with the data within the 
statistical errors. The systematic errors are discussed in the following section. 



7 Systematic Uncertainties 

The procedure used in this paper is tested with approximately eight million hadronic Z° decays 
generated with the JETSET model and passed through a detailed simulation of the OPAL 
detector. The results obtained from the fit, to the reconstructed events, R^ s = 0.224 ± 0.008 
and AFsid, s) = 0.128±0.025, where the errors are statistical only, are found to be in agreement 
with the input values for i?d,s = 0.220 and /li?s(d, s) = 0.100. 

Several sources of systematic uncertainty affect our result. They can be divided into three 
classes: 1) those due to detector effects, 2) those due to model assumptions, and 3) uncertainties 
in the heavy quark contributions. Since these three classes of errors are largely uncorrelated, 
we estimate their individual impact on the branching fractions and asymmetries and add them 
quadratically to obtain the overall systematic error. The individual systematic errors on the 
electroweak observables are summarised in Table 7 for R^ g and in Table 8 for the forward- 
backward asymmetries. The errors are determined by changing each input parameter in turn, 
repeating the analysis, and observing the shifts in the results, as given in the tables. 
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Source of Error 






dE/diX resolution ±5% 


:i:0.0037 


^0.0023 


dE/diX mean ±0.15(7 


±0.0016 


±0.0010 


Kg purity and efficiency ±4.1% 


:i:0.0124 


:f0.0077 


A purity and efficiency ±6.3% 


±0.0003 


±0.0002 


r]f /r]f = 1.00 ±0.02 


< 0.0001 


< 0.0001 


vf^^"^/Vu^, ^^"V^r = 1.00 ±0.02 


^0.0039 


TO.0024 


rif^^'^T]^^ = 1.00 ±0.02 


^0.0022 


TO.0014 


Vu'^'^^/Vd^^^ = 1.00 ±0.02 


TO.0002 


to.oooi 


p(baryon)= p x (1.00 ± 0.02) 


^0.0021 


^0.0013 


= 0.068 ± 0.007 


^0.0005 


^0.0003 


f^^ = 0.101 ± 0.011 


±0.0050 


±0.0031 


/p{p) = 0.088 ± 0.009 


^0.0011 


^0.0007 


= 0.114 ±0.012 


^0.0042 


^0.0026 


fMA) ^ 0.091 ± 0.010 


< 0.0001 


< 0.0001 


0.020 <ec< 0.039 


±0.0015 


±0.0009 


= 0.078 ± 0.004 


TO.0003 


^0.0002 


= 0.039 ± 0.004 


±0.0016 


±0.0010 


/P® = 0.051 ± 0.009 


^0.0011 


^0.0007 


/^s ^ 0.036 ± 0.010 


^0.0033 


^0.0021 


/^(^) = 0.031 ± 0.015 


±0.0002 


±0.0001 


Sub-total Systematic Error 


0.0160 


0.0097 


Rc = 0.158 ±0.010 


< 0.0001 


^0.0037 


i?b = 0.2212 ±0.0019 


< 0.0001 


TO.0008 


Total Systematic Error 


0.0160 


0.0104 


Total Statistical Error 


0.0155 


0.0096 


Total Error 


0.0223 


0.0141 



Table 7: Sources of systematic error and their effects on the measurement of R'^ ^ and R^^g- The 

corresponding errors on R'^ and R^ are twice those of R'^ ^ and -Rd,s and are negatively correlated, 
with the exception of the uncertainties due to R^ and i?b) which are positively correlated. 
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Soiirrp of Flrror 






ySJ-J yJiJu icowlLlLlL'li _l_tJ/0 


zrD DDI S 


3=0 OO'^? 


H / H -T* in pa n — t— n 1 /T 


+n nnriQ 


+0 0095 




-|-L/.UU(J X 


mO 01 ^9 


A purity and efficiency di6.3% 


±0.0002 


±0.0006 
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Table 8: Sources of systematic error and their effects on the measurements of the forward- 
backward asymmetries for E^^ ~ mzo. 
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The following systematic uncertainties are considered: 



Detector sources: 

— hadron identification by dE/dx: As mentioned in Section 3, we apply corrections 
to the apparent charged pion, charged kaon, and proton yields to take into account 
the migration of the true to the apparent particle species and their relative 
efficiencies. The uncertainties in these corrections are estimated by varying the 
widths, a, of the ionisation loss distributions by ±5% and the central values by 
±0.15(7. Since such changes affect all charged particle species similarly, they result 
in smaller uncertainties in the purities of the identified samples than the overall vr^, 
K^, and proton (antiproton) identification efficiencies. 

— efficiencies of Kg and A: The limited knowledge of the Kg and A efficiencies 
and purities affects the discrimination power between up and down type quarks. 
The Kg efficiency is one of the largest sources of systematic uncertainty since the 
difference between Kg and K^*" production separates u and d quark events, since Kg 
tags primarily d and s quarks while K^ tags primarily u and s quarks. 

— other selection criteria: To cross-check for possible biases, we repeat the analysis 
with different maximum polar angles of the tagging particles between 0.6 < | cos6\ < 
0.7, different cuts on the polar angle of the thrust axis 0.7 < | cos6'Thrust| < 0.9, and 
different minimum values of 0.4 < Xp < 0.6. No significant deviations are observed. 
Of the different minimum values of Xp which were studied, a minimum tagging 
hadron momentum Xp > 0.5 was found to result in the smallest combined error. 
The statistical error grows with increasing Xp, while the heavy flavour backgrounds 
(the uncertainty of which is a major source of systematic error) decrease. The 
measurements made using different minimum Xp give consistent results. 

Model uncertainties were addressed in some detail in [^]. Here we summarise the 
important issues: 

— hadronisation symmetries (Eqs. 4—8): The hadronisation relations of Eqs. 4—8 
are stable against variations of fragmentation parameters within 2%. As shown in 
Table 7, the corresponding uncertainties in the branching fractions are small. 

To estimate the uncertainties in the hadronisation relations, we vary several of 
the important fragmentation parameters in JETSETQ. The default values are 
obtained from a fit to event properties such as charged particle multiplicity, the 
Xp distributions for vr^'', K^*", and protons, particle production rates, thrust (and 
other event shape) distributions 0], etc. Uncertainties in the parameters can be 
estimated from the corresponding increase of the of model with respect to 
the data. Requiring a change of at least Ax^=l, we vary the following parameters: 

* the hardness of the fragmentation, achieved by varying the 'b' parameter of the 
LUND symmetric fragmentation function [^| between 0.22 and 0.62, 



* the fraction of strange quarks in the sea, by varying the relevant JETSET 
parameter 7s between 0.27 and 0.36, 



^An alternative hadronisation model is implemented in the HERWIG [ pl[ generator. However, this model is 
in general not SU(2) isospin symmetric due to technical reasons [Q. Therefore, HERWIG is not suited for a 
general study of these hadronisation relations. 
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* the fraction of vector mesons, varied between 0.4 and 0.7 for hadrons with only 
up and down quarks and between 0.2 and 0.6 for those containing strange quarks, 

* the production rates of tensor mesons by varying PARJ(ll) between 0.4 and 0.7, 
PARJ(12) between 0.2 and 0.6, and switching on tensor meson production, 

* the value of Aqcd used in JETSET between 0.25 and 0.29. 

The hadronisation symmetries show no statistically significant variations outside the 
±2% range. 

The relations are also confirmed by a study with the COJETS model which 
invokes independent fragmentation. 

— hemisphere correlations: Dynamic correlations and geometrical effects are taken 
into account by the p parameter, which is found to be independent of the primary 
quark flavour and the tagging hadron type to good accuracy in JETSET. To study 
possible breaking of this universality, we allow p for the proton and the A to vary by 
±2% relative to the overall p, thus allowing for possible differences between baryon 
and meson hadronisation. 

• Heavy flavour contributions: 

— heavy flavour fractions and reliabilities: The estimation of the heavy flavour 
contributions was discussed in Section 4. The largest uncertainties are due to the 
relative amounts of charged and neutral kaons from charm events. The measured 
D° and KgX vs. K^X branching fractions cause the largest systematic 
uncertainties, as can be seen in Tables 7 and 8. It is again the uncertainty in 
the ability to separate up and down quarks with charged and neutral kaons which 
limits the accuracy of the measurements. 

— heavy flavour Standard Model quantities: Variations of the heavy flavour 
branching fractions of the TP and the forward-backward asymmetries are studied and 
found to have only small effects on the measurements of the light flavour electroweak 
parameters (see Tables 7 and 8). Note that the uncertainties on and i?b have a 
negligible effect on the measurements of R'^. 



8 Results and Conclusions 



Taking into account the systematic uncertainties, the final results (assuming R'^ = R'^) are: 

R'^, = 0.371 ±0.016 (stat.) ±0.016 (syst.) and (29) 
K = 0.258 ±0.031 (stat.) ±0.032 (syst.) . (30) 

The Standard Model predictions are 0.359 and 0.282, respectively. The i?^ have the advantage 
of being largely independent of Rc and i?b- If we assume the world average values of Rc and 
i?b from [^, we find: 

i?d,s = 0.230 ±0.010 (stat.) ±0.010 (syst.) and (31) 
= 0.160 ±0.019 (stat.) ±0.019 (syst.) . (32) 
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For the forward-backward asymmetries, we determine: 



AFB{d,s) = 0.068 ± 0.035 (stat.) ± 0.011 (syst.) and (33) 
ApBin) = 0.040 ±0.067 (stat.) ±0.028 (syst.) (34) 

for centre-of-mass energies near the Z° peak. The measurement of AFB{d,s) is consistent with 
the result obtained by the DELPHI collaboration |^ using a different method which is more 
dependent on hadronisation models. Note that the measurements R'^ ^ and R'^ are completely 
negatively correlated and those of ^^^^(d, s) and Apsi^) more than 90% positively correlated. 
All of the measurements are in agreement with the predictions of the Standard Model. 

We interpret these results in terms of right and left handed Z°gg couplings g'^^ and g'j^^ of 
down and strange quarks. Here g1 = Pq{Il — tq sin^ Q^s) and g\ = —pqCg sin^ 9es, Pq and sin^ ^eff 
being the effective electroweak parameters, I3 the third component of the weak isospin, and 
the quark electric charge in units of the magnitude of the electron charge. The left and right 
handed couplings are obtained by solving the equations 

R' (4")^ + ial'r .35^ 

A^B(d,s)(mzo) = (36) 

4 {9r r + {9l r 

with the electron coupling Ae = 0.1466 ± 0.0033 0. The right and left handed couplings 
of up quarks are assumed to be fully anticorrelated to the down/strange ones such that the 
denominator of Eq. 35 is equal to the Standard Model expectation. In addition, the energy 



dependence of the off-peak asymmetries to lowest order was parametrised according to : 



Giis) = QlQ) I r ±2 QeQfVeVf Re (x^Xzo) + {v^ + al) {vj + a}) \ Xzo |' (38) 
G3{s) = 2 QeQf aeQf Re {x*^xzo) + 4:VeaeVfaf \ Xzo \^ (39) 

with 

Xzo = ^— ^- (40) 

s — m^o + iTzos/niio 

and 

X'y = ^ , H^ = -0.0593 ± 0.0007. (41) 

The axial and vector couplings are related to the left and right handed couplings through 
aj = gi ^ 9r and vj = gi + g{i- The values of the axial and vector couplings for the electron 
are taken from p[. 

The measurements are corrected for initial photon radiation and QCD effects using the 
JETSET model including initial state photon radiation. Both hard gluon and photon 
emission are suppressed due to the requirement of a high-Xp particle. The corrections are 
negligible for the measurements of Rq and for the asymmetries amount to ±0.002 due to QCD 
effects and ±0.004, ±0.002 and ±0.017 due to initial state photon emission at centre-of-mass 



22 



energies of 89.6, 91.2 and 92.9 GeV, respectively, which are small compared to the statistical 
error. 

We find g'^^ = — 0.441o;q9 and g'^^ = 0.13l;o;i7 in agreement with the Standard Model values 
of —0.424 and +0.077, respectively. The results are shown in Fig. 5 together with various 
confidence level regions. The fit yields a of 2.6 for two degrees of freedom. The couplings enter 
quadratically into the on-peak observables but linearly in the 7-Z° interference contribution 
for the off-peak asymmetries. Therefore, only the latter measurements allow discrimination 
between the possible signs of gL and gji. The solution for a negative g^ leads to an acceptable 
X^=3.0 and cannot be excluded. Solutions with a positive gL are disfavoured and can be 
excluded with 76% and 83% confidence levels for negative and positive g^, respectively. 

Whereas results on down quarks exist from lepton-nucleon scattering, the present analysis 
is the most direct measurement of the strange quark couplings. Departing from the assumption 
of the equality of the couplings of down and strange quarks, we fix R'^ and determine R'^, 
yielding the results presented in Fig. 6. Fixing R'^ to the Standard Model value results in 
i?s = 0.392 ± 0.043 (stat.) ± 0.045 (syst.). The dependence of R'^ on R'^ is linear and can be 
parametrised as dR'^/dR'^ = —1.83. Note that the value of R'^ is completely anticorrelated 
to -Rg. Similar results can be obtained for Afb{S') and ApBivi.) for various fixed values of 
^FB(d) as shown in Fig. 7. For Apsid) fixed to the Standard Model value of 0.100, we find 
Apb{s) = 0.075±0.028 (stat.)±0.008 (syst.) and ^^.^(u) = 0.086±0.030 (stat.)±0.021 (syst.))%. 
In this case, ^^^(s) and Apsi^) are only +31% correlated. Again, the dependence of the 
results on ApBid) is linear and can be parametrised by dAFB{s)/dAFBid) = +0.32 and 
dApBiyi) / dApB^d) = +1.42. The statistical and systematic errors on Aj?s(u) are independent 
of Apsid), while those of Apb{s) scale linearly with the value of Apb{s) such that the relative 
error is constant. 

A measurement of all of the individual light-fiavour electroweak parameters would require 
that the number of unknown fiavour tagging efficiencies, r]^, be reduced. This might be possible 
in the future by using additional hadronisation symmetries which, however, may introduce 
a greater model dependence. 
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Figure 1: Two particle separation power {dE/dxi — dE / dx2) / {cr{dE / dx)) in units of mean 
dE/dx resolution as a function of particle momentum for tracks with | cos^| < 0.7 and at least 
20 measured samples, obtained from data. The curves given are for pion/proton, pion/kaon 
and kaon/proton separation. 
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Figure 2: Invariant mass tt+tt" of Kg candidates with Xp > 0.5 in the data (points with error 
bars) compared to Monte Carlo (histogram) which is normahsed to the same number of Z° 
events. Note that the details of the total Kg production rate in the Monte Carlo are not 
relevant to this analysis. The arrows indicate the mass range in which Kg candidates are 
selected as high-Xp tags. The different mass resolutions in the data and Monte Carlo are taken 
into account in the efficiency determination. 
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Figure 3: Invariant mass p7r~ of A candidates with Xp > 0.5 in the data (points with error bars) 
compared to Monte Carlo (histogram) which has been normahsed to the same number of Z° 
events. The arrows indicate the mass range in which A candidates are selected as high-Xp tags. 
Note that the A fragmentation function is too hard in JETSET, resulting in an overestimation 
of the high-Xp production of A baryons, and that the details of the total A production rate in 
the Monte Carlo are not relevant to this analysis. 
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Figure 4: ApBidjS) and Apsi'^) versus the centre-of-mass energy, -Ecmj where the errors 
are statistical plus systematic. The statistical errors for different centre-of-mass energies 
are uncorrelated. The systematic errors are correlated between centre-of-mass points. The 
Standard Model predictions for mtop = 175 GeV and mmggs — 300 GeV are shown as the sohd 
curve for AFBid^s) and the dashed curve for Apb{u). 
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Figure 5: Left vs. right handed coupUngs of strange quarks, assuming identical values for down 
quarks. The 68% and 95% confidence level regions are also shown. The points correspond to 
the measurement and the Standard Model expectation. 
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Figure 6: R'^ = Rg/ (Rd+Ru+Rs) determined for various fixed values of R'^ — R^/ (Rd+Ru+Rs)- 
The result obtained by enforcing the equality of R'^ and R'^ is given as the star, and the Standard 
Model value is shown as the solid point. The hatched region shows the total error on R'^ for a 
given R'^. Note that R'^ = l- R'^- R'^. 
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Figure 7: Apsi^) and Apsi'^) determined for various fixed values of ^^^^(d). The results 
obtained by enforcing the equality of the down-type light flavour couplings are shown as the 
triangle and the star and the Standard Model values are shown as the square and circle for 
u and s quarks, respectively. The single and double hatched regions show the total error on 
^fb(u) and Apb{s), respectively, for a given Apsid). 
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